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ABSTRACT 

We present the first maps of the surface magnetic fields of a pre-main sequence binary 
system. Spectropolarimetric observations of the young, 18 Myr, HD 155555 (V824 Ara, G5IV 
+ KOIV) system were obtained at the Anglo- Australian Telescope in 2004 and 2007. Both 
datasets are analysed using a new binary Zeeman Doppler imaging (ZDI) code. This allows 
us to simultaneously model the contribution of each component to the observed circularly 
polarised spectra. Stellar brightness maps are also produced for HD 155555 and compared to 
previous Doppler images. 

Our radial magnetic maps reveal a complex surface magnetic topology with mixed polar- 
ities at all latitudes. We find rings of azimuthal field on both stars, most of which are found to 
be non-axisymmetric with the stellar rotational axis. We also examine the field strength and 
the relative fraction of magnetic energy stored in the radial and azimuthal field components 
at both epochs. A marked weakening of the field strength of the secondary star is observed 
between the 2004 and 2007 epochs. This is accompanied by an apparent shift in the location 
of magnetic energy from the azimuthal to radial field. We suggest that this could be indicative 
of a magnetic activity cycle. We use the radial magnetic maps to extrapolate the coronal field 
(by assuming a potential field) for each star individually - at present ignoring any possible 
interaction. The secondary star is found to exhibit an extreme tilt (w75°) of its large scale 
magnetic field to that of its rotation axis for both epochs. The field complexity that is apparent 
in the surface maps persists out to a significant fraction of the binary separation. Any interac- 
tion between the fields of the two stars is therefore likely to be complex also. Modelling this 
would require a full binary field extrapolation. 

Key words: Stars: pre-main sequence - Stars: magnetic fields - (Stars:) binaries: spectro- 
scopic - Polarization - Stars: imaging - Stars: coronae - 



1 INTRODUCTION 

In the past decade our understanding of the surface magnetic 
topologies of stars beyond our Sun has greatly increased. This 
is mostly due to the applicatio n of the Zeeman-Doppler imaging 
(ZDI) technique ( ISemellll989l) which allows us not only to map 
magnetic regions on stellar surfaces, but also to determine the ori- 
entation of the field lines in these regions. From combining obser- 
vations in circularly polarised light with the positional information 
(using Doppler imaging principles), we can determine whether a 
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particular region on the stellar surface is composed of mainly ra- 
dial or azimuthal field. ZDI has now been applied to many rapidly 
rotati ng active stars with some surprising results (see lDonati et al.l 
1 20071 for a recent review). 

The magnetic maps of active G and K dwarfs show regions 
of strong surface azimuthal field. For many of these targets the re- 
covered azimuthal field strengt hs are in ex cess of the radial field 
(e.g. for AB Dor and LQ Hya, D onati et al j 2003), thus suggesting 
that the magnetic energy is concentrated in the toroidal field rather 
than the polodial field. This is quite different from what we find on 
the Sun where most of the resolved magnetic structures are radi- 
ally orientated field lines. Current solar dynamo theory tells us that 
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strong azimuthal field should only be found at the shear layer be- 
tween the radiative core and convective envelopes (the tachocline). 
lDonatietalJ ill999h therefore suggest that the observation of strong 
surface azimuthal field on these rapidly rotating stars leads to the 
conclusion that the Q, dynamo component is not confined to the 
tachocline but instead may be distributed throughout the convec- 
tive envelope. 

The only tidally locked binary star to date for which published 
magnetic maps are available is the evolved primary star of the RS 
CVn system HR 1099. A strong latitudinal dependence on the po- 
larity of the radial field maps is found and HR 1099 often exhibits 
a unipolar cap (Donati et al. 2003). Also present are strong, and of- 
ten com plete, axisy mmetric rings of azimuthal field on the stellar 
surface. iPetit et al.l ( 12004^ also confirmed the existence of an anti- 
correlation between the polarity of the radial and azimuthal fields. 
This is not observed on young rapidly rotating stars but is present 
on the Sun. Due to the evolved nature of RS CVn primary com- 
ponents, it is difficult to interpret their magnetic maps with respect 
to those of young rapidly rotating single stars. Therefore in an at- 
tempt to establish the relative importance of tidal locking and youth 
on surface magnetic field distributions, we present observations of 
the young pre-main sequence binary system HD 155555. 

The HD 155555 (V824 Ara) system is composed of a G5 IV 
primary and a KO IV secondary and has an orbital period of 1.68 
d. This makes it a good target for observing over a five day ob- 
serving run, as outlined i n j{2l HD 155555 was first discovered as 
a close binary system by iBennett et alj 1 1967h. It w as later clas- 
sified as an RS CVn binary jStrassmeier et alj|l993h based upon 
the activity and binarity of the system. Pasguini et alj jl99lh sug- 
gested that HD 155555 was more likely a young pre-main sequence 
binary. This was based upon the high Li (6708 A) abundance (re- 
visited here in i]3.2t and the pre sence of a very active dMe compan- 
ion (LDS587B). More recenflv lStrassmeier & Ricd ( l2000h derived 
an age of 18 Myr using the Hipparcos distan ce of 31 pc and the 
pre-m ain-sequence evolutionary tracks of D' Antona & Mazzitellil 
(Il997h . 

HD 155555 is a particul arly active system wit h an X-ray lu- 
minosity of 2.7xl0^''erg dDempsev et alj[l993h . both compo- 
nents show Call H & K core emission and filled in Ha. Therefore 
given its proximity and activity it is unsurprising that HD 155555 
has been the subject of a number of recen t papers. In 199 6 it wa s 
the target of a multi-wavelength study by iDempsev et alj ( 1200 ih . 
The authors detected considerable flare activity, with several flares 
being observed simultaneously in the UV and radio wavelengths. 
As part of this campaign, optical spectroscopy was also obtained 
and sur face brightness images wer e produced using Doppler imag- 
ing by IStrassmeier & Rica ( I2OO0I) . An earlier Doppler image of 
H D 155555 produced from data obtained in 1990 was published 
by iHatzes & KtirsteJ ill 999*1. In addition, a magnetic field detec- 
tion has also been r eported for both components of HD 155555 by 
lDonatietalJ ( ll997h . 

In this paper we present new surface brightness maps (Sj3) and 
the first magnetic maps of HD 155555 ( i]4.1b . To achieve this we 
have developed a binary ZDI code in order to model the contri- 
bution of each star to the combined Stokes V profiles that are ob- 
served during conjunction phases. We outline the details of this in 
section f|4] The more recent, and numerous, 2007 observations of 
HD 155555 are described first. Then in SQwe present a smaller set 
of observations take in 2004. In obtaining magnetic maps for both 
components of a binary system we open up many interesting pos- 
sibilities to explore the effects of binarity on the magnetic fields 
of both stars. Possible interaction between the stellar fields are of 



interest because it will determine many of the systems X-ray prop- 
erties and the location of stellar winds. In ^we use the magnetic 
maps recovered of both the primary and secondary components in 
order to perform an initial analysis of the likelihood of interaction 
between the two stars. We discuss our results in ij7]and present our 
conclusions in 38] 



2 OBSERVATIONS AND DATA REDUCTION 

Spectropolarimetric observations of HD 155555 were made at the 
Anglo-Australian Telescope (AAT) using the University College 
London Echelle Spectrograph (UCLES) which is fibre fed by the 
SemelPol visiting polarimeter ( ISemeletalJll993l) mounted at the 
Cassegrain focus. In total eleven uninterrupted nights of observa- 
tions were secured on our target star HD 155555 in 2007 from 
March 30 to April 09. In this paper we only use data collected on 
five of these nights, from March 31 to April 04. The remaining 
six nights of observations are used in a forthcoming publication 
to study the surface rotation properties of both stars. The 1.68 d or- 
bital period means that good phase coverage can be obtained in five 
nights, corresponding very closely to three stellar orbits. 

The instrumental set-up is the same as that used on many 
previous Zeeman Do ppler imagi ng runs and is described in detail 
by previous authors ( Semel et al . 1993, Donati et al. 2003). Using 
SemelPol both left- and right-hand circular polarised light are cap- 
tured simultaneously on the EEV2 CCD detector. Using the 31 gr 
mm~^ grating full wavelength coverage between 4377 and 6819 A 
is achieved. A resolution 70,000 (i.e. 4.3 km s~^) is obtained us- 
ing this set-up. In order to supress all spurious polarisation signals 
(down to first order) due to instrumental drifts or temporal varia- 
tions in stellar spectra (e.g. resulting from the orbital motions of 
the components in a binary system), circular polarisation (Stokes 
V) observations are made from a sequence of four exposures. Be- 
tween each exposure the quarter-waveplate in SemelPol is rotated 
by +1- 45°, causing a switch of the circular polarisation carried in 
each fibre. The four sub-exposures are then combined using op- 
timal extraction techniques as implemented by the Echelle Spec- 
tra Reductio n: an In teractive Tool (ESpRIT) software developed by 
Donati et ail h997l) . In the process the normal calibrations pertain- 
ing to spectral data are performed and the spectra are wavelength 
calibrated. As HD 155555 is such a bright target (V = 6.9 mag) we 
were able to use 200 s exposures which meant that, after readout 
time, the total time taken for each sequence of four exposures (and 
so each Stokes I and V spectrum) was approximately 15 minutes. 

In total 53 spectra were obtained over the five nights and are 
fairly evenly spread over the orbital phase of the binary. The signal- 
to-noise (S/N) ratio of both our Stokes I and V spectra in the peak 
order vary with the weather conditions between 70 and 230 with 
an average of around 160 per 1.9 km s~^ pixel. In order to greatly 
boost the signal-to-noise the technique of Least-Squares D econ- 
volution (LSD) is used as described by iDonati et alj jl99 7) and 
implemented in ESpRIT. This effectively combines the absorption 
line profiles of many thousands of photospheric lines, resulting in 
a single high S/N ratio line profile. For HD 15555 5 we use a line- 
list for a KO star produc ed by the VALD database jPiskunov et al.l 
ll995llKuDkaetalJ ll999'). All strong chromospheric lines (such as 
the Balmer series) are excluded. Approximately 2800 lines are used 
for LSD which result in a S/N ranging from 2500 to 10,500 for 
Stokes V spectra and between 680 and 1200 for Stokes I spectra. 
We take advantage of the strong telluric lines present in the spectra 
and the LSD technique in order to correct for drifts in our wave- 
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Table 1. Orbital and physical parameters for the HD 155555 system. 



Element 


Unit 


Value 


Source 


Orbital: 








P 


(days) 


1.6816463 


adopted^ 


To 


(HJD) 


2446997.9102 


adopted^ 


<I>0 




0.752474 


DoTS 


e 




0.0 


adopted 


7 


(kms-i) 


3.72±0.02 


DoTS 


? 




0.935±0.001 


DoTS 


Ki 


(kms-l) 


86.4±0.1 


DoTS 


K2 


(kms-l) 


94.7±0.2 


(Ki and q) 


i 


(degs) 


50 (52 adopted) 


DoTS 


{v sin i)i 


(kms-i) 


34.9 


DoTS 


{v sini)2 


(kms-l) 


31.3 


DoTS 


Physical (assuming i = 52°): 






mi 


(Mo) 


1.054 


DoTS 


m2 


(Mq) 


0.986 


DoTS 


Ri 


(R@) 


1.47 


DoTS 


R2 


(Rq) 


1.32 


DoTS 


TcB,l 


(K) 


5300±100 


SME 


Teff,2 


(K) 


5050±100 


SME 


log 91 




4.05±0.1 


SME 


log 92 




4.10±0.1 


SME 



^ from^teassmeier_&Rice QOOO). 



length calibration over the course of the nig ht. lDonatietai]i2003l) 
show that this calibration results in radial velocities which are sta- 
ble to better than 0.1 km s~^. 



3 SYSTEM PARAMETERS AND BIUGHTNESS MAPS 

We use the Doppler imaging code 'DoTS ' jCollierCameronlll997h 
to map the surface brightness distribution of the surfaces of both 
stars using the Stokes I intensity spect ra. This code use s maxi- 
mum entropy, as implemented by .Skilling & BrvanI l ll984 'l. to re- 
turn the simplest possible surface map at a given level of x^. A 
two- temperatu re model is used to describe the surface of each star 
jCoU ier-Camer on & Unru h 1994). Spectroscopic look up tables are 
created to model the contribution from each element on the stel- 
lar surface. Slowly rotating reference stars are fitted with Gaus- 
sian profiles which we use to model both the local intensity profiles 
of the photosphere and the spot, for each star. It has been shown 
bylUm'uh & Collier Cameron ( 1995) that using synthetic Gaussian 
profiles instead of those from slowly rotating stars has little effect 
on the resulting brightness maps. While the obvious advantage of 
using synthetic Gaussians is that they are noiseless. For the pri- 
mary star we use temperatures of 5300 K and 4050 K for the pho- 
tosphere and spot respectively. Similarly the cooler secondary star 
is modelled using a photospheric temperature of 5050 K and a spot 
temperature of 3787 K. These temperatures come from the spectro- 
scopic analysis described in ^3.2\ 

3.1 Determining system orbital parameters 

In addition to producing surface brightness maps we can use DoTS 
to determine the fundamental parameters of the HD 155555 system. 
This is achieved by finding the values of each system parameter that 
give the best fit to the stellar profiles (i.e. that mini mise x^) after 
a fixe d number of maximum entropy iterations (see iBames et al.l 
I2OOOI) . In a binary system the additional orbital parameters greatly 



increase the computer time required to achieve a global solution. 
We approached this problem by carefully selecting combinations of 
two parameters to be minimised for simultaneously. In practice we 
used a pool of 30 computers to simultaneously run different com- 
binations of the variables. Several iterations were made through the 
pairs of parameters until a global convergence on a solution with 
= 0.6 was achieved. At this point further changes in the parameters 
were of the order of their uncertainties as estimated from the curva- 
ture of the surfaces. Values for system parameters are collated 

in Table[T] 

We use the ephemeris of IStrassmeier & Ricd bOO Ol) which is 
based on over twenty years of radial velocity data. We attempted 
to derive our own value for the period of the system using our 
2007 dataset alone and fo und that it was consistent with that of 
IStrassmeier & Ric3 j2000l) . however with naturally larger uncer- 
tainty. The value of 0.752474 we obtain for the phase offset ($0) 
reflects the traditional eclipsing binary definition that we use in 
DoTS. Here phase zero occurs at stellar conjunction with the pri- 
mary star furthest away from the observer (i.e. normally primary 
eclipse in an eclipsing binary system). We make the simplifying 
assumption that both components of HD 155555 have rotational 
axes perpendicular to the orbital plane. Unfortunately determining 
the inclination of the binary orbit for non-eclipsing systems is as 
problematic as the rotation axis of single stars because the Doppler 
imaging process is only weakly sensitive to this parameter. We de- 
rive a value of 50° (using the minimisation technique described 
above) which is consistent with the value of 52° obta ined by both 
IStrassmeier & Ric3 ( l200(]f) and IPasguini et alj ( Il99ll) . For consis- 
tency with these works we therefore also adopt a value of 52° . 

3.2 Atmospheric and evolutionary parameters 

A subset of our highest signal-to-noise spectra are selected to deter- 
mine the atmospheric parameters using spectral synt hesis. We ap- 
ply th e technique of tomographic separation (see Bagn uolo & GiesI 
I1991I) to disentangle the two components using the orbital pa- 
rameters derived in the last section. We then calculated and fitted 
synthetic spectra to each compo nent with the help of the Spec- 
troscopy Made Easy (SME) code jValenti & Piskunovlll996t) . This 
code uses least-squares minimization to determine the atmospheric 
parameters that best describe the observed spectra. For th ese cal- 
culatio ns we used stellar model atmospheres from the iKuruczl 
1 1993|) grid, and atomic line data obt ained from the VALD database 
('Pisk unov eTZIl 19951; fKupka et alj [ 1999). In our calculations we 
assumed solar abundances, a microtubulence of 1.75 km s~^ and a 
macrotubulence of 2.0 km s^^, all typical values for a G5~K0 class 
pre-main-sequence object. 

For our analysis we selected two wavelength regions, one 
around the temperature and gravity sensitive Na I D lines, and one 
between 6000 and 6200 A, a region that contains a large number 
of metal lines. The line ratios in the latter regions provide indepen- 
dent leverage on the effective temperatures. We find the effective 
temperatures and gravities of the two stars to be very similar (pri- 
mary: Teff,! = 5300 ± 100 K and log gi = 4.05 ± 0.1, secondary 
Tcff,2 = 5050 ± 100 K and log 32 = 4.10 ± 0.1). The full set of 
parameters is listed in Table [T] 

Within their uncertai nties, the values we find are consistent 
with the values quoted in Pasguini et al. I( ll99lh . The main differ- 
ence we find is that the surface gravity of the s econdary star is sig- 
nificantly lower than the value of 4.5 quoted bv lStrassmeier & Ricd 
{2OOO) and consequently much closer to the value obtained for the 
primary star. This is also expected from the orbital parameters of 
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the system, since g ~ GM/R^ and thus: 

A(log5) = log 31 - logp2 = log[A^i/M2 * {R2/Rif] 

Substituting M1/M2 = K2/ K\ and assuming synchroneous rota- 
tion, R2/R1 — (f sin 1)2/(1; sini)i this gives: 

A(log5r) = \og{K2/Ki * [(usini)2/(«sin j)i]^) 

Using the values from Table[T] one finds A(logg) = —0.05. 

We also re-examine the lithium abundance to provide an in- 
dication of youth. Using the atmospheric parameters above, we 
find logn(Li) = 3.16 (primary) and 3.15 (secondary). Our 
analysis takes account of the non-LTE corrections published by 
ICarlsson et all ( 1 19941) that were not available to jPasguini et al.l 
( 199 ih . These lithium abundances are slightly below the primor- 
dial values of 3.2-3.3, a nd suggest a lithium age of 0~30 Myr 
dSestito & RandiclJ |2003) . This is in good agreement with the 
age of 18 Myr derived by [Strassmeier & Rice (2000) from the 
pre-main sequence tracks of IP'Antona & Mazzitellil ( Il997h . We 
therefore agree with the c onclusions of Pasguini et al.l 1 199lh and 
IStrassmeier & Ricel ( l2000h that HD 155555 is a young pre-main 
sequence binary system and not an evolved RS CVn object. 



3.3 Brightness images 

Now that the best parameter set for HD 155555 has been found, 
we can attempt to model the deviations (bumps) in the line pro- 
file due to the presence of stellar spots. Unlike earlier in this sec- 
tion where we were interested in the parameter set that minimises 
X^, when generating final spot maps a target value is set such 
that the Maximum Entropy iterations converge before we start over- 
fitting the data and producing artificially noisy maps. In Fig.ll2lwe 
show the 53 observed Stokes I LSD profiles and the correspond- 
ing fits we obtain using DoTS that correspond to a reduced x^= 
0.7. A value for below unity is often achieved when perform- 
ing Doppler imaging using intensity spectra (Stokes I), it simply 
indicates that the S/N values from the LSD process are slightly un- 
derestimated. The profiles have been ordered by phase so the binary 
motion can clearly be seen and individual spots are easily traced. 
From the LSD profiles alone we can see that the flat bottomed shape 
of the rotation profiles of both stars indicates the presence of polar 
spots. 

In Fig. |2] the brightness maps produced by the Doppler imag- 
ing process are shown. We choose to leave these maps labelled with 
phase, as opposed to the more conventional stellar longitude, so that 
the locations of spots can be more easily compared to the binary or- 
bital phase. Note that for the primary star phase (j) — 0.0 faces the 
secondary star, while for the secondary star phase (fi = 0.5 faces 
the primary. The first obvious characteristic of both images are the 
large polar spots. On the primary star individual darker spots can 
be made out at latitudes between 60° and the pole. A band of low- 
latitude spots between (10 - 40°) is also present and covers three 
quarters of the stellar circumference. While we have been able to 
resolve some structure in this band it is likely that in reality it is 
composed of smaller, unresolved spots. Latitudes between 40 - 60° 
are remarkably devoid of spot coverage producing the bi-modal 
plot of mean spot filling factor versus stellar latitude shown in the 
right panel of Fig. [2] 

The secondary star has a large offset polar spot that is centred 
at phase (j> = 0.4 and reaches down to a latitude of 60°. A large 
spot group at (/) = 0.2 is the lowest latitude feature reaching down 
to 20° and is clearly seen in the LSD profiles of Fig. [121 A weak 
pair of mid-latitude spots is also recovered at phases (p = 0.7 and 



(j) = 0.8. In Fig.[3]we show the brightness maps projected on to the 
stellar surfaces so that the locations of spots can be compared to the 
binary orbital phase. We discuss the starspot distributions further in 
f|7]where we also compare them to previously published maps. 



4 ZEEMAN-DOPPLER IMAGING OF BINARY SYSTEMS 

All previous magnetic maps have been produced using ZDI codes 
designed to work on single stars. Therefore the few maps of the 
primary stars of binary systems that have been produced to date 
were obtained by performing the initial step of correcting the LSD 
profiles for the orbital motion. The data can then be treated as that 
of a single star. This is a common procedure for modelling of both 
Stokes I and V data by codes that are not capable of modelling 
the contribution from both stars. Due to the continuum present in 
Stokes I spectra this is particularly difficult. The binary conjunction 
phases are either not observed, or the part of the primary rotational 
profile that overlaps with the secondary star is masked out, or the 
secondary profile is removed by subtraction of a Gaussian. The sit- 
uation when using Stokes V data is considerably easier however, 
this is because the continuum is unpolarised. For RS CVn binaries 
the magnetic field strength of the secondary star is so small (due in 
part to the ratio of the stellar luminosities) that its contribution to 
the combined Stokes V signatures can often be assumed negligible. 
It is this assumption that is used to produce the many magnetic im- 
ages of the prim ary star in the RS CVn binary system HR1099 (e.g. 
IPetit et al.ll2004h . 

In the case of the near equal mass ratio binary HD 155555, 
the stars are so similar that the above options are not viable. For 
approximately one third of the orbital period the rotation pro- 
files of the stars overlap to some extent, as illustrated in Fig. [T2] 
If one chose to simply not observe the conjunction phases then 
there would be two phase gaps each of duration A(j) ~ 0.16. 
These would at the very least reduce the level of detail recovered 
in the images and may even produce imaging artifacts. Gaps in 
the phase coverage is worse for magnetic than brightness imag- 
ing due to the fact that we ascertain the nature of the field (e.g. 
radial or azimuthal) by the behaviour of the Stokes V signature as 
it moves through the line profile. As we shall see, the magnetic field 
strengths on the two stars are comparable. Therefore the only effi- 
cient way to proceed is to develop a code capable of performing 
full ZDI of a binary star system. 

We base our new ZDI code on that of lHussain et~ZI ilOOdt 
which itself was originally spawned from the DoTS code, we there- 
fore name the new code 'ZDoTS'. One of the main reasons why 
the binary version was not previously implemented was due to the 
computational demands that this would place. The memory require- 
ments of ZDoTS exceeds that of the DoTS by a factor of three due 
to the necessity of storing each of the magnetic vectors (the ra- 
dial, azimuthal and meridional fields). ZDoTS uses array compres- 
sion and takes advantage of the general increase in recent computer 
memory in order to overcome this problem. 

When the two binary components are in conjunction we model 
the result as a linear superposition of the two Stokes V signatures 
from both stars. Just as in Doppler imaging observations, such 
phases alone do not provide information about the relative contri- 
bution to the observed Stokes V spectrum from each star. However, 
when combined with observations where the two components are 
separated in velocity space then knowing how the two Stokes V 
signatures combine helps greatly to constrain the geometry of the 
magnetic field maps. 
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Figure 1. The resulting 2007 LSD profiles are plotted along with their uncertainties. The solid line is the maximum entropy regularised fit we obtain to the 
data from Doppler imaging. Careful examination of the Une profiles reveals the presence of spots, seen as raised bumps in the line profiles. Profiles are ordered 
by orbital phase (listed to the right of each profile) and clearly show the motion of both binary components. For reference, the primary star is the profile seen 
generally at negative velocities for phases = 0.0 — 0.5. Included in the on-line version of this paper is a plot showing the starspots in more detail. 



As in the single star ZDI code of iHussain et alj ( llOOOl) . we 
model the local Stokes V line profiles using a Gaussian profile 
(which is assumed constant over the stellar surface) and we as- 
sume that the weak field approximation holds. Three consecutive 
Stokes V spectra of HD 155555, shown in Fig.|4l illustrate partic- 
ularly well what can happen when the stars are in conjunction. The 
Stokes V signatures from the two stars can be seen moving closer 
together between the top and bottom spectra of Fig. |4] Then dur- 
ing the conjunction phase in the bottom spectrum we observe, by 
chance, the exact superposition of the troughs of the two profiles 
which produces a large amplitude combined signature. As Fig. |4] 
shows, the ZDoTS code has successfully modelled this. We note 
that the three spectra shown in Fig. |4] were not obtained in the five 
nights considered in this publication but from the subsequent five 
nights of data (specifically, April 05). The ZDoTS code was tested 
by creating synthetic Stokes V spectra from an input map with spots 
located at different 

In order to test ZDoTS further we used the surface magnetic 
maps of the single, KO d warf, AB D or (u sini = 91 km s ~^) from 
1996 and 1999 (Donati & Collier Ca meron 1997 and Do nati et al.l 
l2003h . We 'painted' one of these maps on to the primary star and 
the other onto the secondary and assumed all the same orbital and 
physical parameters of the HD 155555 system (from Table[Tll. Syn- 
thetic spectra were then generated, using the forward module of 



ZDoTS, of comparable signal-to-noise and at each of the same 
53 observational phases as our real data. This resulted in Stokes 
V spectra of the binary system including the tangled conjunction 
phases. The spectra were then used as input into ZDoTS and we 
attempted to reconstruct the magnetic maps. The results were com- 
pared with the perfor mance of the original (single star) ZDI code of 
iHussain et al.l jlOOOl) at recovering the magnetic maps of each star 
in turn. The results showed that both codes recovered essentially 
the same maps. When the images were collapsed in latitude the 
maximum difference between the two magnetic fluxes was always 
less than 5%. Importantly, there was no visible phase (longitude) 
dependence on the accuracy of the reconstructed maps. This illus- 
trates that ZDoTS is as capable of reproducing surface maps for 
binary systems as for single stars, at least for well phase sampled 
data. To examin e further tests of ZDoTS using synthetic data see 
iDunstond |2003) . 



4.1 Magnetic maps of HD 155555 

In Fig.|5]we show the 53 Stokes V LSD profiles along with the max- 
imum entropy regularised fits corresponding to x'^= 1-0. As in Fig. 
I12lwe order the Stokes V profiles by phase and again the binary or- 
bital motion can be clearly seen. Unlike in Fig. [121 it is easy in Fig. 
[5]to spot the few spectra which have substantially poorer signal-to- 
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Figure 2. Doppler images of the surface brightness distribution of both components in 2007, showing the location of the prominent stai' spots. The tick marks 
at the bottom of each main panel mark the phases of observation. To the right of each map is a graph showing the mean spot filling factor as a function of 
stellar latitude. Note that for the primary star this is significantly bimodal. 



noise. These occur around phases (j) = 0.52 and = 0.77 and still 
provide some constraint to modelling the surface magnetic field. 

The strengths and complexity of both the primary and sec- 
ondary Stokes V signatures varies with orbital phase. From a vi- 
sual inspection of Fig. [5] it appears from the Stokes V profiles that 
the primary star often possesses a more complex field with either 
more changes of polarity and/or different field orientations than the 
secondary star. This is illustrated by the greater number of sign re- 
versals seen throughout the line profile of the primary star (as many 
as six seen around phase (j) — 0.15). The Stokes V signature from 
the secondary star all but disappears around phase (j> = 0.6. 

The reconstructed magnetic images are shown in Fig. (6) We 
show the radial, azimuthal and meridional field components. The 
density and coverage of the observations ensures that there are no 
imaging artifacts that can result from poor phase sampling. The 
maximum amplitude of both the radial and azimuthal field is ap- 
proximately twice as large on the primary as on the secondary star. 
Note that ZDI is only sensitive to high latitude meridional flux on 
stars with moderate to high inclination angles. At low latitudes the 
meridi onal map suffers considerable cross-talk with the radial map 
jPonati & Brow ni l 1997b . A comparison of the recovered low lati- 
tude meridional field (bottom panels Fig. [6]l with the radial field 
(top panels Fig.[6]l shows this to be the case. 

The radial field map of the primary star shows many small 



regions of flux of alternating polarities. The level of surface detail 
recovered is in fact near the limit that our observations of a star 
with a relatively small v sini (~ 35km s^^) could produce. As an 
interesting comparison, we note that the radial field maps of the 
primary star are similar in the level of detail they recover to those 
produced of AB Dor (actual v sini=91km s~^) once spun-down to 
V sin i=35km s~^, as was performed during the testing of ZDoTS 
described earlier. 

In Fig.|7]we show the radial magnetic field painted on to the 
surface of each star. This illustrates the field structure near the polar 
regions better than the mercator projection (Fig. |6j. It also allows 
us to consider the distribution of flux on the star with respect to 
the binary orbital motion. We discuss the magnetic maps further in 
i|7.2| where we will also compare them to those obtained in 2004. 



5 2004 OBSERVATIONS 

A small number of observations were also acquired at the AAT of 
HD 155555 during 2004 September 23 - 28. The instrumental set- 
up and data reduction procedures were identical to those described 
in ij2] The exposure time of each individual exposure was 300 s 
(compared to the 200 s of the 2007 observations). All orbital param- 
eters of HD 155555 were re-calculated for the 2004 observations 
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Figure 3. The 2007 Doppler brightness images, in Fig. |2] ai'e shown projected onto the surfaces of each star at the four quadrature phases. The unspotted 
photosphere of the secondary star appears darker to the primary star due to the difference in stellar temperatures. 
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Figure 4. Three example Stokes V LSD profiles are shown to illustrate the 
ability of our code to reproduce the complex supeiposition of Stokes V 
signatures that arise when the stars are in conjunction. Note the two Stokes 
V minima (one from each star) that move closer together due to orbital 
motion between the top and middle spectra. They then combine in the third 
spectrum to produce the largest amplitude signal we observe during our 
campaign. We note that the signal-to-noise of the top two spectra shown 
here are below average. 



using the DoTS code and the minimisation technique (see t]3.1t . 
The mass ratio (g) and semi-major velocity amphtude (Ki) were 
found to be consistent with the 2007 values (within the 2004 values 
uncertainties), therefore we adopted the higher precision 2007 val- 
ues shown in Table [T] The phase offset (0o) and the radial velocity 
(7) were also re-calculated for the 2004 dataset. A phase offset of 
00 = 0.752606 (cf. 00 = 0.752474, Table [T) and a radial veloc- 
ity of 7 = 3.77 km s"^ (cf. 7 = 3.72 ± 0. 02 kms~\ Table[B 
were f ound. We note that the value found bv IStrassmeier & Ricd 
1I2OOO') for 1996 was considerably larger, 7 = 5.9±0.2kms"\ 
One possible source for changes in the apparent radial velocity of 
the HD 155555 is the dMe companion, LDS587B. However, a ~2 
km s^ ^ change is probably too large to be attributed to the well sep- 
arated LDS587B. Future observations will be required to determine 
if long-term trends exist in both 0o and 7. 

Over the six nights a total of seven Stokes V spectra were 
taken. Two of these were taken close to each other on the same 
night and on the sixth night a very similar phase as was observed 
on the first night was obtained. Therefore together the Stokes V 
spectra effectively only sample five unique phases. The largest gap 
in the phase coverage is A(j) ~ 0.25 between phases = 0.51 - 
0.76. The signal-to-noise in the peak order ranged from 180 - 220 
on six of the nights and was 290 on September 25. In Fig. [8] we 
show the 2004 Stokes V spectra and the fits. From these alone it 
is possible to ascertain that the field structure has changed on both 
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Figure 5. The 2007 resulting observed Stokes V LSD profiles are plotted along with their uncertainties. The solid line is the maximum entropy regularised fit 
we obtain to the data from Zeeman Doppler imaging. Profiles are ordered by orbital phase (listed to the right of each profile) and clearly show the motion of 
both binary components. 



stars. The secondary star generally shows a greater peak-to-peak 
amplitude that it did in 2007 and is now comparable to the primary 
star. 

The relatively small number of observations in 2004 compared 
with the 2007 dataset result in lower quality maps with less detail 
and increased phase smearing. This is especially true for the bright- 
ness maps shown in Fig. [9] Here the weak low latitude spots are 
only recovered by pushing the maximum entropy regularised fits to 
the point where noise starts to appear. Despite this fact the 2004 
data can still be compared with the 2007 data to provide valuable 
insight into the evolution over the intervening 2.5 year period, this 
is discussed in i|7.1| 

The magnetic maps we recover are also shown in Fig. |9] The 
meridional magnetic maps show considerable cross-talk with the 
radial map. This much more severe in 2004 than the 2007 maps 
(Fig.|6]l. T his is to be expected due to the poor phase coverage at this 
epoch (see lDonati & Brownll 19971 for further details). We note that 
on close comparison of the 2004 maps with those in Fig. |6] (2007) 
it is apparent that regions in the 2004 maps are more stretched out 
in latitude than those in 2007. We suspect that this is a result of the 
poor phase sampling of the 2004 data. In order to assess the impact 
this has on the 2004 magnetic maps, we created a subset of seven 
spectra from the 2007 dataset that were nearest in phase to those 
of 2004. We then proceeded to use only these spectra in the ZDl 
process to re-produce the 2007 maps. The results of this procedure 
were that essentially the same maps were recovered however with 
increased latitudinal smearing present. This also led to a reduction 
by nearly 50% in the field strengths of the recovered magnetic re- 
gions. 



These results are not surprising given that the Doppler imag- 
ing process is inherently more sensitive to longitude than latitude 
position and it is only by repeatedly observing the same features as 
they move through the stellar rotational profile that we can accu- 
rately determine their latitude position. Also, as mentioned earlier, 
ZDI is only able to ascertain the orientation of the field by observ- 
ing how the Stokes V signature behaves as it rotates into and out 
of our line of sight. Therefore the poorer the phase sampling the 
greater the cross-talk between the radial, azimuthal and meridional 
fields. Interestingly though, because this affects both stars equally, 
the ratio of the amplitudes of the field strengths recovered between 
the primary and secondary stars remained similar for this subset of 
2007 data as for the whole dataset. We postpone further discussion 
of the 2004 magnetic maps to Wl ■2\ 



6 CORONAL FIELD EXTRAPOLATIONS 

The surface radial magnetic map of each star can be extrapolated 
to provide us with an initial idea of the coronal field struc t ure us - 
ing the coronal X-ray modelling technique of Jardin e et al.l (|2003), 
based around a code originally developed bv ivan Ballegooiien et al.l 
( Il998l) . At present we treat each star separately, making no attempt 
to model the interaction between the two stars that may occur. 
While this obviously places limitations on the conclusions we will 
be able to draw from the field extrapolations, it will have only min- 
imal effect on either the small scale structure near to the stellar sur- 
face or the global orientation of the magnetic axis. This technique 
uses the "Potential Field Source Surface" method which assumes 
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Figure 6. The 2007 magnetic surface maps produced from the Zeeman Doppler imaging process. Maps of the primary star are shown on the left, while those of 
the secondary are on the right. The top panels show the radial magnetic field, the middle shows the azimuthal maps and the bottom panels show the meridional 
field. In all maps the neutral grey colour represents zero magnetic field, white shows positive polarity and black shows regions of negative polarity. Note that 
the greyscale contrast range is different between the primary and secondary maps, as illustrated in the key to the immediate right of each map. 



two boundary conditions, the first being that the radial field is equal 
to that recovered on the radial field map at the stellar surface and 
secondly that the field becomes purely radial at the source surface. 
The implication of this is that the magnetic field is forced open by 
the outward pressure of coronal gas at the source surface radius. 

The further we get from the surface of a star the simpler the 
field geometry becomes. Eventually the outermost groups of closed 
field lines represent the lowest order, dipolar, field. From these we 
can establish whether or not the magnetic axis is aligned with the 
stellar rotational axis. The largely unsee n lower hemisphere con- 
tains little magnetic flux. As discussed bv lJardin^ et al.l h002l) . the 
nature of the field in the unseen hemisphere can influence the coro- 
nal geometry of the low latitude field although the effect is much 
less apparent for the high-latitude field. The presence of a binary 
companion may also change the connectivity of the largest-scale 
coronal field lines, since these may possibly connect the two stars. 
The orientation of the magnetic axis of the field, however, is not 
affected by either of these factors since it is determined purely by 



the surface field of the observed hemisphere. In order to illustrate 
the coronal topology we choose the source surface to be at a large 
enough radius for each star that the higher-order components of the 
field have died away and only the simplest, dipolar, structure re- 
mains. This is at a radius of 3.4 R, for the primary and 6.8 R, for 
the secondary component. The true locations of the source surface, 
and so the radii of the last closed field lines, will of course depend 
upon the density of X-ray emitting coronal gas. The value obtained 
fo r HP 155555 of l og (rie) = 10.7 ± 0.23 [cm^^] determined 
bv lNess etaLl ( l2004l) would suggest that the true source surface is 
probably sm aller, perhaps compar able to the values determined for 
AB Dor by iHussain et al. I ( l2007i) . We note however that on sin- 
gle rapidly rotating stars stellar prominences have b een observed at 
severa l stellar radii above the surface. For example. iDunstone et al.l 
( l2006h found that Speedy Mic had a complex prominence system 
located at twice the height above the stellar surface as the calculated 
co-rotation radius. This leads to the requirement of large closed 
magnetic loops in order to enforce co-rotation. In comparison the 
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Figure 7. The 2007 radial magnetic maps, in Fig. |6] are shown projected onto the surfaces of each star at the four quadrature phases. Note that we use a 
greyscale contrast range of ±400 G which is a compromise of the different scales used for each star in Fig.|6] 



co-rotation radius of HD 155555 is in effect the orbital separation 
of 7.5 -R0 (by virtue of the fact that this is a tidally locked system). 

In Fig.[TO]we show a representative sample of field lines from 
both regions close to the surface of the star and the larger scale field. 
As we described in the last section the radial surface maps of both 
stars, but the primary star in particular, show a complex magnetic 
field with many small magnetic regions. We can thus expect there to 
be significant structure at relatively small heights from the surface 
of each star. Such regions will produce small closed loops. This 
indeed is what we see from the coronal field extrapolations in Fig. 
[Tol 



7 DISCUSSION 

7.1 Distribution of stellar spots 

In addition to the September 2004 and April 2007 brightness maps 
presented in this work, t here are two previo u sly pu blished Doppler 
images of HD 155555. Hatzes & Kiirster (1999) present obser- 
vations taken in September 1990 and Strassmeier & Rice (200(?) 
those obtained in May 1996. These four well spaced epochs pro- 
vide a valuable opportunity to study the evolution and stability of 
the spot groups in this tidally locked binary system. The details 
of the different Doppler imaging codes, the phase coverage and 



signal-to-noise of observations mean that we have to be careful in 
our interpretation of the strength and sizes of individual features 
that we see. However a general comparison of the distribution of 
spots should be possible. 

Conside ring the primary star firs t, all four maps show a polar 
spot. In 1990 Hatzes & Kurstej (Il999) observed a significantly de- 
centred polar spot that was tilted in the direction of the secondary 
star The other three epochs show a more centrally located polar 
cap that extends down to a latitude of -1-60° . The well sampled 2007 
dataset shown in Figs.[2]and[3]suggest that the polar cap is not ho- 
mogeneous but is composed of a number of large dark spots. All 
four epochs show a lack of mid-latitude spots and a band of low 
latitude spots at approximately -1-30°. This bi-modality in the lati- 
tudinal distribution of spots seems to be a consistent feature of the 
G5IV primary. This is reminiscent of the bri ghtness map s found 
for G dwarfs in the young a Persei cluster bv lBames et all lll998l) 
and more recently of the G dwarf HD 171488 by Marsden et al.l 
( 2006). All these stars show spots emerging near the polar or at low 
latitudes but with a distinct lack of spots at intermediate latitudes. 
Theoretical flux-tube modelling by Granzer et al. ( 2000) predict a 
bimodal spot distribution for very young pre-main sequence stars. 
IStrassmeier & Ric3 ( l2000h considered these models for the particu- 
lar case of HD 155555 and found that the predictions did not match 
the observed spot patterns. Instead, meridional flows may explain 
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Figure 8. The observed 2004 Stokes V profiles from the LSD procedure 
are plotted along with their uncertainties. The solid line is the maximum 
entropy regularised fit we obtain to the data from Zeeman Doppler imaging. 
Profiles are ordered by orbital phase (listed to the right of each profile). 

how flux is transported to the poles. Or else two distinct flux emer- 
gence mechanisms may be in operation. 

Hatzes & Kiirster ( 1999) revealed a string of such low-latitude 
spots on the anti-facing hemisphere (i.e. the side of the primary 
facing away from the sec ondary star). This was also reported by 
IStrassmeier & Ric3 ( |2000|) . although their maps recover spots cov- 
ering a larger range of stellar longitudes and therefore also revealed 
a number of spots on the trailing hemisphere (the one that always 
follows behind with respect to the binary orbital motion). Both our 
maps for 2004 and 2007 reveal a band of low-latitude spots. The 
more reliable 2007 maps show that this covers three quarters of the 
stellar circumference and shows a preference for spots to be located 
on the facing and trailing hemispheres. The gap that breaks what 
would otherwise be a complete ring of low-latitude spots occurs be- 
tween phases <j> = 0.2 - 0.5. This corresponds to the leading/facing 
hemisphere of the primary star and is thus in total contrast to the 
1990 map. We are therefore left to conclude that there are no appar- 
ent stable 'active longitudes' on the primary star of the HD 155555 
system. 

With the exception of 2004, the secondary star does not have a 
true polar cap but instead ha s a hig h-latitude spot that just touches 
the pole. iHatzes & KtirsteJ ( 1 19991) found such a spot tilted to- 
wards the primary star (a t phase (j> ~ 0.55) in 1990, whereas 
IStrassmeier & Ric3 j2000l) recovered a similar feature but on the 
trailing hemisphere (at phase <j} ~ 0.30) of the secondary star in 
1996. Our 2004 spot map shows the secondary star to possess a 
true (centrally located) polar spot. In 2007 we again find a high- 
latitude spot that this time is centred on phase c/> ~ 0.4 and so is 
again tilted in the direction of the primary star. Given our relatively 
poor understanding of how polar spots form, it is possible that the 
shifting location of the polar spot may be linked to the extreme tilt 
of the magnetic axis found in Sj6] 

All epochs reveal spots at mid and low-latitudes on the sec- 
ondary star. In particular a large low-latitude (30°) spot group is 
seen at all epochs between phases <^ = 0.1 — 0.3 (trailing hemi- 



sphere). As this appears to be a consistent feature of the secondary 
star maps it may indeed indicate an active spot longitude. 

7.2 Surface magnetic topology and dynamo processes 

The magnetic maps we produced in the last sections (Figs. |6] and 
^ are the first of a pre-main sequence binary system. Both compo- 
nents of HD 155555 possess large regions of strong azimuthal field 
near to the stellar surface. These are of equivalent strength (in fact, 
often stronger than) the recovered radial field. Strong azimuthal 
field is recovered on all active cool stars that have been mapped 
so far using ZDI. This is in contrast to the relatively small amounts 
of horizontal field found on the comparatively inactive Sun. So- 
lar azimuthal field is thought to be buried deep in the base of the 
convective e nvelope where dyna mo activity is genera ted. Previous 
authors (e.g. IPonati et al.ll2003l and IPetit et"^l2004l) have there- 
fore postulated that the presence of strong azimuthal field, at the 
stellar surface of these rapidly rotating stars, suggests that their dy- 
namo may be active throughout the convection zone. While this 
suggestion is contrary to the classical understanding of the gen- 
eration of large-scale fields in slowly rotating (solar-type) stars, 
it has recent l y gain ed support from the numerical simulations of 
iBrown et alj J2007l) . These authors find that global-scale toroidal 
and poloidal magnetic fields can be built and maintained in the 
bulk of the convection zone of rapidly rotati ng stars, despite the 
presence of turbulent convection. In particular lBrown et al.l i2007h 
find strong ordered torodial fields which may explain the observa- 
tion, as in HD 155555, of strong surface azimuthal field. We can 
now compare the maps of both components of HD 155555 to each 
other and look for evolution between the 2004 and 2007 epochs. 



7.2.1 Primary star 

The radial field map of the primary star shows many small regions 
of flux of alternating polarities. A visual inspection of the 2007 pri- 
mary radial maps reveals no real latitude dependence, with both 
positive and negative field at all latitudes and no preference at the 
pole. In 2004 the maps do show a slight preference for positive 
polarities, especially at mid to high latitudes (60 - 80°). The radial 
field maps are integrated over all longitudes in Fig.[TT]to test for the 
presence of a latitudinal polarity dominance. In doing so consider- 
able cancellation of the flux occurs due to the already described 
complex nature of the maps. However, Fig. [TT] shows that we are 
left with a preference for positive radial flux at all latitudes and for 
both epochs, with a peak flux at 55°. 

The azimuthal field maps of the primary star reveal a consider- 
ably simpler field structure than the radial maps. This is particularly 
apparent in 2007 where the surface topology can essentially be de- 
scribed by just two regions. A large region of positive field is seen 
encircling the star at high latitudes, centred on (~75°). This points 
down to lower latitudes (~60°) at phase (j) = 0.8. It therefore re- 
sembles a ring of positive flux that is tilted by 15° with respect to 
the rotational axis. Due to the fact that the ring crosses over the ro- 
tational pole, part of it then appears as a negative polarity region. 
The other region is of negative field at lower latitudes. In 2004 the 
azimuthal maps appear more complex that in 2007, with more mix- 
ing of positive and negative polarities. The general latitude distri- 
bution is similar to 2007, with positive field at high latitudes and 
negative field at lower latitudes. However, the field appears more 
axisymmetric in 2004 with both rings of flux being better centred 
on the rotation pole. Again, in Fig.[TT]we integrate the azimuthal 
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Figure 9. Brightness (top) and magnetic maps (middle - radial field, bottom - azimuthal field) are shown for September 2004 observations. Note that unlike 
2007 maps, the greyscale contrast range of the magnetic maps is the same between the primary and secondary maps, as illustrated in the key to the immediate 
right of each map. 



maps over all longitudes to find that the primary star indeed shows 
a dominant negative polarity at all latitudes below 55°, but then 
switches to positive field at high latitudes. The similarity of the az- 
imuthal flux distribution in both 2004 and 2007 suggests that this 
is a relatively stable feature of the primary star. 

The meridional field maps show two very high latitude re- 
gions of opposite polarity in both 2007 and 2004. This is exactly 



what you would expect to result from a homogeneously oriented 
ring of horizontal field that crosses the pole. The ZDI process will 
recover part of the field as azimuthal (as discussed above) and part 
as meridional. The low latitude features (as already outlined in i]4.1l 
and SjSjl are the result of cross-talk from the radial field (particularly 
severe in 2004). Therefore, the azimuthal and high-latitude merid- 
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Figure 10. Coronal field extrapolations are shown for the primary star (left panels) and secondary star (right panels) for 2007 data (top panels) and 2004 data 
(bottom panels). These particular phases were chosen as they best illustrate the large scale field orientation. The primary is shown at phase </< = 0.8 (and so 
near to the trailing hemisphere), while the secondary is at phase <j> = 0.9 (nearest to the anti-facing hemisphere). Open and closed field lines are represented 
by blue and white lines respectively. The radial magnetic maps are shown on the stellar surfaces with blue corresponding to negative, green to zero and red to 
positive magnetic flux. 

field at phases (f} ~ 0.6 — 0.9. However, we see a positive ring of 
field at high latitudes, again tilted (like the azimuthal field of the 
primary star in 2007) away from the rotational axis by ~ 15° and 
towards phase (j) = 0.3. Fig.[TT]shows that negative azimuthal field 
dominates at latitudes below 45° then the field becomes positive at 
higher latitudes. Similarly in 2004, a ring of non-axisymmetric pos- 
itive azimuthal field can be seen at high-latitudes. Fig. [TT] shows a 
preference for positive azimuthal field at all latitudes. Both epochs 
have a strong, low latitude, region of negative azimuthal field cen- 
tred on phase ~ 0.5 (the hemisphere facing the primary star). 
The meridional maps of the secondary star show two high-latitude 
regions of opposite polarity. They are therefore similar to those of 
the primary star and we thus refer the reader to the discussion of 
the primary star's meridional maps in the last section. 



ional flux appears to be attributable to the toroidal field, while the 
radial and low-latitude meridional flux is part of the poloidal field. 



7.2.2 Secondary star 

The secondary star has a considerably simpler radial field map than 
the primary star. Like the corresponding map of the primary star, 
the radial field map appears to have no latitude dependence on the 
polarity of the field. Indeed when we collapse the maps in longitude 
in Fig.[TT]we achieve virtually total cancellation of the positive and 
negative fluxes. In 2007 the residual trend is for positive polarity 
at latitudes lower than 55°and negative flux at higher latitudes. In 
2004 there is a slight indication of preference for positive polarity 
for both field components. 

The maps of azimuthal field of the secondary star are unusu- 
ally complex. In 2007 we recover low latitude flux of mixed polar- 
ity, with a complex region composed of both negative and positive 
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Figure 11. The latitudinal flux distribution is shown for each star on both the 
2004 and 2007 epochs. In each panel the solid line represents the quadratic 
magnetic flux, the dashed line shows the azimuthal field and the dash-dot 
line depicts the radial field. The value of a given quantity is the integration 
over all longitudes and has been corrected for the relative pixel area at that 
particular latitude - hence the reduction in all fluxes observed at the pole. 



7.3 Magnetic regions, polar spots and axis misalignment 

To summarise the last section; we find evidence for rings of az- 
imuthal field which are non-axisymmetric but yet show a clear con- 
stant polarity structure. The trend for both stars of HD 155555 is to 
have a high-latitude ring of positive azimuthal field which is tilted 
with respect to the rotational axis and then for this to be under- 
lined by negative field at low latitudes. Interestingly the rings of 
positive azimu thal field closely follo w the polar spots discussed in 
fTTTj Donati ( 199^ and lDonati et al . ( 2003) found a similar pattern 
between the azimuthal field and polar spots on HR 1099. 

The similarity of the 2004 and 2007 coronal extrapolations 
shown in Fig.[TO]is striking. On both epochs the secondary star ex- 
hibits an extreme tilt of its magnetic axis. From plotting the neutral 
polarity line at the source surface we calculate that the magnetic 
axis is approximately 75° misaligned with the rotational axis. The 
primary stars magnetic field is more closely aligned to that of its 
rotational axis. We find a ~30° misalignment to that of the rota- 
tional axis in 2004 and a larger ~45° misalignment in 2007. This 
may be related to the fact that, as discussed in the last section, the 
ring of surface azimuthal field shows a tilt by approximately 15° 
for the primary star in 2007 but more of an axisymmetric ring in 
2004. So it appears that both the radial and azimuthal fields were 
both more axisymmetric in 2004. It therefore seems likely that the 



surface azimuthal field is linked to the orientation of the large scale 
field. 

7.4 Relative field strengths and magnetic energy 

In the 2007 maps shown in Figs.|6]and|9]the magnetic field strength 
scale needed to be different for each star and yet for 2004 they could 
be the same. The average quadratic field intensity can be calculated 
for both stars at each epoch. This can also be examined graphically, 
as a function of latitude, from the solid line in Fig. [TT] We find 
for the well sampled 2007 epoch an average field strength (inte- 
grated field modulus over the whole stellar surface) of 68 G for the 
primary star and 43 G for the secondary star. So in 2007 the pri- 
mary star appeared to have a field strength 50 % greater than that 
of the secondary star. For the 2004 maps we find a field strength of 
52 G for the primary star and 61 G for the secondary. It therefore 
appears that the field strength of the secondary star was consider- 
ably stronger in 2004, both in absolute terms (42% stronger than in 
2007) and relative to that of the primary star. 

As mentioned in ^the comparatively poor phase sampling of 
2004 data must be taken into account. We found that using only a 
subset of the 2007 data (choosing the same number of spectra avail- 
able in 2004) resulted in an almost 50% reduction in the recovered 
flux. In addition the observed flux was spread out in stellar latitude. 
The effect of this can be seen in Fig. [TT] where the total quadratic 
flux (solid line) shows a broader distribution in the 2004 panels 
than those of 2007. The effect of fewer phases of observation is 
that we effectively miss (do not recover) flux from the stellar sur- 
face. Therefore the 16 G weaker field strength of the primary star 
in 2004 is likely due, at least in part, to this effect. However, re- 
duced phase sampling certainly cannot explain why the secondary 
star has a field strength 18 G (42%) stronger in 2004 than in 2007. 
In fact it suggests quite to the contrary, the real difference in the 
flux between 2004 and 2007 must have been even greater. 

Another quantity which we can calculate in an attempt to un- 
derstand the above observations is that of the magnetic energy. This 
can be calculated for each of the radial and azimuthal fields. We can 
then quote the energy stored in each field component as a fraction 
of the total energy. The meridional field is found to contain be- 
tween 8 - 12 % of the total magnetic energy (depending upon the 
star and epoch). It is therefore small compared to the radial and az- 
imuthal contributions and we do not discuss it further We find that 
in 2007 the azimuthal field of the primary star stores 58% of the to- 
tal energy with 34% in the radial component. This is similar to the 
case of AB Dor and LQ Hya where the azimuthal field is found to 
preferentially store more of the field energy than the radial field at 
most epochs iD onati et alj2003l) . In the case of the secondary star, 
in 2007 the radial field dominates by storing 48% of the magnetic 
energy while the azimuthal field has 43%. 

When we perform the same calculation for the 2004 maps we 
find that the azimuthal and radial fields store equal energy (44%) 
on the primary star, corresponding to a decrease in the relative az- 
imuthal field strength. However on the secondary star the opposite 
appears to have occurred, the azimuthal field is found to contain 
54% of the energy and the radial field with 38%. We again note 
that the poor sampling of the 2004 dataset may have an impact on 
these results. Observations with very poor phase sampling can be 
prone to cross-talk between the radial and azimuthal field compo- 
nents. This would lead to a reduction in the contrast between the 
two field components and may partly explain the observation of 
more equal percentages for the primary star. Also if, by chance, a 
particularly strong region of field (be it azimuthal or radial) was 
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missed by our observations then this would be reflected in the bal- 
ance of the magnetic energy. The situation would be considerably 
worse in the case of poor phase coverage. If the observations were 
very poorly spaced then there would be little constraint on the flux 
for large parts of the stellar surface. While the 2004 observations 
are poorly sampled they are at least reasonably well spaced. We 
therefore conclude that the relative strengths of the field compo- 
nents for 2004 should be treated with caution and as an indication 
only. 

In summary, it appears that the magnetic field of the secondary 
star has evolved considerably between 2004 and 2007. The aver- 
age field strength has become much weaker (both in absolute value 
and relative to the primary star) and the azimuthal field component 
appears to have weakened relative to that of the radial. This is ac- 
companied by an apparent change in the azimuthal topology from a 
comparatively simple field (characterised by a ring of positive flux) 
in 2004 to a more c omplex distribution in 2007. If, as suggested by 
iDonati et alj ( l2003h . the strong azimuthal field observed on active 
stars is closely related to the toroidal component of the large-scale 
dynamo field then the observed behaviour may represent the decay 
phase of an activity cycle. There is insufficient evidence to make a 
similar comment about the primary star at this time. 

7.5 Comparing HD 155555 with other stars 

The HD 155555 system is particularly interesting as it is a young 
binary system and is therefore a potentially very promising target 
for disentangling the effects of youth and binarity. As a close binary 
system, the components of HD 155555 are tidally locked just like 
the primary stars in evolved RS CVn binaries, such as HR 1099. 
However, with an age of around 18 Myr (see ^3.2\ the system is 
still a pre-main sequence object and therefore younger than the well 
studied zero-age main-sequence ra pidly rotating sing le stars, AB 
Dor and LQ Hya (50 - 100 Myr, e.g. ljanson et alj2007r) . This fact is 
reflected in the bloated size of both components of HD 155555. De- 
spite the fact that AB Dor and LQ Hya are of approximately equal 
mass to HD 155555 we find, in iJHI] radii of 1.5 and 1.3 Rq for 
the primary and s econdary star respe ctively (c.f. « 1 Rq for AB 
Dor and LQ Hya, ' Donati et al.ll2003 f). This will lead to increased 
convection zone depths of both components of HD 155555 relative 
to AB Dor and LQ Hya. Therefore, somewhat ironically, both the 
evolved HR 1099 binary (R=3.7 Rq) and HD 155555 have convec- 
tion zone depths greater than AB Dor and LQ Hya. Consequently 
if convection zone depth uniquely determines the observed mag- 
netic surface topologies, then our attempt to compare the effects of 
binarity and youth will be ultimately inconclusive. 

As already mentioned, the observation of strong surface az- 
imuthal field is common to HD 155555, AB Dor, LQ Hya and HR 
1099. Beyond this, HR 1099 is characterised by easily identified 
rings of azimuthal field that are normally axisymmetric. A similar 
ring of azimuthal field is often observed on AB Dor but only at 
high-latitudes. We find that both components of HD 155555 show 
rings of azimuthal field, however they are often non-axisymmetric. 
HR 1099 also possesses a latitudinal dependence on the polarity 
of the radial field with a unipolar polar cap. In contrast the young 
rapid rotators AB Dor and LQ Hya show mixed polarities of ra- 
dial field at all latitudes, including the pole. At first glance this is 
essentially what is seen for HD 155555 but when the latitudinal 
distribution of flux is plotted a clear trend is seen. It can be seen in 
Fig.[TT]that the primary star (in both 2004 and 2007) has a bias to- 
wards positive radial flux at all latitudes but increasing so towards 
high-latitudes. IPetit et alj l2004h confirmed the presence of an anti- 



correlation of the latitudinal distribution of azimuthal and radial 
field on HR 1099. This is not seen on AB Dor or LQ Hya but is 
present on the Sun. Fig. [TT] shows tentative evidence of a similar 
anti-correlation for HD 155555 in 2007 for the secondary star and 
for the primary star at latitudes below 55° in both 2007 and 2004. 
However, this does not hold for the secondary star in 2004 where 
both the radial and azimuthal field is positively orientated. Further 
observations of HD 155555 will be needed to establish long-term 
trends. 



7.6 The coronal fields of HD 155555 and binary interaction 

In Sj6]we used the recovered radial field maps to extrapolate the 
coronal structure of both stars. Due to the small binary separation 
of 7.5 Rq (~ 5.5 i?*, see Table[B there is no doubt that the mag- 
netospheres of the two stars will interact. The question is: what will 
the nature of the interaction be? Due to the relative complexity of 
the surface radial field maps of both stars, the extrapolated coronal 
field is also correspondingly complex. As Fig. [TO] shows, at loca- 
tions between the stars where the two fields would interact, both 
fields retain sufficient complexity that their interaction will not be 
simply modelled by that of two dipolar fields. 

It is clear that the next step will be to perform a full binary ex- 
trapolation, where the radial field maps of both stars are considered 
simultaneously. This will model which field lines will preferentially 
connect with the other star rather than to regions on the same star. 
The upper limit on the longevity of such 'binary field lines' will 
largely be defined by the relative differential rotation rates of each 
star. The surface rotation properties of HD 155555 will be the sub- 
ject of a forthcoming paper. 



8 CONCLUSIONS 

We have produced the first simultaneous magnetic maps of both 
components of a binary system. These show that the technique of 
Zeeman Doppler imaging can be adapted to the binary case. This 
opens up the potential to study many short period pre-main se- 
quence and main sequence binary systems where both components 
have Stokes V detections. The HD 155555 system provides an in- 
teresting first insight into the magnetic fields of a young binary 
system. In general, we find that the resulting magnetic maps ap- 
pear to share features in common with both young rapidly rotating 
stars and the evolved primary stars of RS CVn binaries. The most 
surprising result is the severe misalignment of the secondary stars 
magnetic and rotational axes. The impact of this can be seen on the 
radial and azimuthal maps. Also interesting is the significant weak- 
ening of the secondary star's magnetic field strength from 2004 to 
2007. Especially as this appears to be associated with a decline in 
the proportion of magnetic energy stored in the azimuthal field. 

Once a binary potential field extrapolation code has been com- 
pleted it will be possible not only to more reliably model the coro- 
nal structure of both stars but also their interaction. This will allow 
us to predict the locations of stellar winds and determine which 
field lines are capable of containing hot X-ray emitting gas. We can 
therefore look forward to learning more about the mass loss rates 
from binary systems and modelling the X-ray properties of bina- 
ries, including flare loop lengths and energies. 
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Figure 12. ONLINE ONLY: The 2007 Stokes I LSD profiles and fits are shown after subtraction of the model immaculate (unspotted) stellar profiles. The 
solid line is the maximum entropy regularised fit we obtain to the data from Doppler imaging. Profiles are ordered by orbital phase (listed to the right of each 
profile). For reference, the primary star is the profile seen generally at negative velocities for phases <f) = 0.0 — 0.5. The residuals in the centre of the stellar 
profiles are starspots and are mostly well modelled. We note ho wever that poor agree ment is found at the extremities of the rotational profiles of both stars. 
This is similar to that observed for the contact binary AE Phe bv lBames et al.l j2004l) . Here the authors attributed the mis-match to inadequate modelling of 
centre-to-limb brightness variation. A more sophisticated treatment of limb darkening beyond the linear limb darkening used here is required. 



